Astronomy & Astrophysics manuscript no. arxivversion 


©ESO 2012 


February 22, 2012 





The Spitzer Spectroscopic Survey of S-type Stars 

K. Smolders 1 *, P. Neyskens 2 '**, J.A.D.L. Blommaert 1 , S. Hony 3 , H. Van Winckel 1 , L. Decin 1 ***, S. Van Eck 2 , 
G. C. Sloan 4 , J. Cami 5 , S. Uttenthaler 1 - 6 , P. Degroote 1 , D. Barry 4 , M. Feast 7 ' 8 , M.A.T. Groenewegen 9 , M. Matsuura 10 , 
J. Menzies 8 , R. Sahai 11 , J. Th. van Loon 12 , A.A. Zijlstra 13 , B. Acke 1 '***, S. Bloemen 1 , N. Cox 1 , P. de Cat 9 , 
M. Desmet 1 , K. Exter 1 , D. Ladjal 1 , R. 0stensen ! , S. Saesen 114 , F. van Wyk 8 , T. Verhoelst 1 -***, and W. Zima 1 

1 Instituut voor Sterrenkunde (IvS), Katholieke Universiteit Leuven, Celestijnenlaan 200 D, B-3001 Leuven, Belgium 
e-mail: kristof . smolders@ster.kuleuven.be 

2 Institut d'Astronomie et d'Astrophysique (IAA), Universite Libre de Bruxelles, C.P.226, Boulevard du Triomphe, B-1050 
Bruxelles, Belgium 

3 Service d'Astrophysique, CEA Saclay, 91 191 Gif-sur-Yvette, France 
e-mail: sacha . hony@cea . f r 

4 Astronomy Department, Cornell University, Ithaca, NY 14853-6801, USA 

5 Department of Physics and Astronomy, University of Western Ontario, London, Ontario N6A 3K7, Canada 

6 Department of Astronomy, University of Vienna, TiirkenschanzstraBe 17, 1 180 Vienna, Austria 

7 Astrophysics, Cosmology and Gravity Centre, Astronomy dept., University of Cape Town, 7701 South Africa 

8 South African Astronomical Observatory, PO Box 9, Observatory, 7935, South Africa 
rj^ 9 Koninklijke Sterrenwacht van Belgie, Ringlaan 3, B-l 180 Brussel, Belgium 

10 Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK 

11 Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91 109, USA 
04 12 Lennard- Jones Laboratories, Keele University, Staffordshire ST5 5BG, UK 

13 Jodrell Bank Centre for Astrophysics, The University of Manchester, School of Physics & Astronomy, Manchester M13 9PL, UK 
' j ' 14 Observatoire de Geneve, Universite de Geneve, Chemin des Maillettes 51, 1290 Sauverny, Switzerland 

Ph 

J/} February 22, 2012 

abstract 

Qh 

Context. S-type AGB stars are thought to be in the transitional phase between M-type and C-type AGB stars. Because the composition 
^ of the circumstellar environment reflects the photospheric abundances, one may expect a strong influence of the stellar C/O ratio on 

the molecular chemistry and the mineralogy of the circumstellar dust. 

Aims. In this paper, we present a large sample of 87 intrinsic galactic S-type AGB stars, observed at infrared wavelengths with the 
i i Spitzer Space Telescope, and supplemented with ground-based optical data. 

Methods. On the one hand, we derive the stellar parameters from the optical spectroscopy and photometry, using a grid of model 
atmospheres. On the other, we decompose the infrared spectra to quantify the flux-contributions from the different dust species. 
^ Finally, we compare the independently determined stellar parameters and dust properties. 

Results. For the stars without significant dust emission features, we detect a strict relation between the presence of SiS absorption in 
the Spitzer spectra and the C/O ratio of the stellar atmosphere. These absorption bands can thus be used as an additional diagnostic 
for the C/O ratio. For stars with significant dust emission, we define three distinct groups, based on the relative contribution of certain 
dust species to the infrared flux. We find a strong link between group-membership and C/O ratio. Furthermore, we show that these 
groups can be explained by assuming that the dust-condensation can be cut short before silicates are produced, while the remaining 
free atoms and molecules can then be used to form the observed magnesium sulfides or the carriers of the unidentified 13 fim and 
20 /im features. Finally, we present the detection of emission features attributed to molecules and dust characteristic to C-type stars, 
such as molecular SiS, hydrocarbons and magnesium sulfide grains. We show that we often detect magnesium sulfides together with 
molecular SiS and we propose that it is formed by a reaction of SiS molecules with Mg. 
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1. Introduction between M- and C-type stars in which the C/O ratio of stars 

changes from the solar value to ratios larger than unityP] In the 
The dredge-up of carbon in asymptotic giant branch (AGB) stars spectra i classification scheme, S-type stars are defined by the 
plays a key role m the chemical evolutio n from oxygen-rich M- preS ence of absorption bands of ZrO and LaO molecules in the 

optical spectra. The transition from oxygen-rich to carbon-rich 



type stars to carbon-rich C-type stars (Iben & Renzini 1983 



|Herwig||2005| |Habing & Olofsson||2004) . S-type stars are of- marks an i mportant chemical transition in the AGB evolution 



ten considered to be an intermediate phase of AGB evolution Following the chemical pathways relevant in the circumstel- 

lar environment, CO is the first molecule to form, locking up all 
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of the available carbon or oxygen, whichever is less abundant. 
Subsequently, in M-type stars (C/O < 1) almost all C atoms will 
be consumed by CO and oxygen-rich molecules such as TiO, 
SiO, and H2O will be present. In carbon-rich stars (C/O > 1) 
the dominant molecules are CO, CH, C2, and C2H2. However, 
it is important to realize that this is a simplified picture of the 
circumstellar chemistry and only valid under the assumption of 
local thermal equilibrium (LTE). Due to non-equilibrium chem- 
istry effects, oxygen-rich stars could display molecules that are 
typical for carbon-rich stars and vice versa ( Cherchneff|2006 1 



The chemical and physical properties of circumstellar dust 
grains are linked to the chemical composition of the stellar atmo- 
sphere. In an oxygen-rich environment, silicates and oxides are 
formed, while graphite, amorphous carbon, and silicon carbide 
are formed in carbon-rich environments. When the C/O ratio is 
very close to unity, both O and C are almost completely con- 
sumed in CO, in which case the dust condensation sequence is 
not well studied nor understood. Theoretical models predict the 
presence of FeSi and metallic Fe, but this still needs observa- 
tional confirmation (Ferrarotti & Gail 2002). Furthermore, some 
molecules and dust species that contain neither C nor O are also 
limited to certain chemical environments, for example magne- 
sium sulfide dust and SiS molecules are typical for carbon-rich 
environments. This indicates that the C/O ratio does not only in- 
fluence the carbon- and oxygen-abundance, but indirectly effects 
the entire chemical network, of which only a couple of molecular 
products are detected. 

The chemistry of the circumstellar material is key to un- 
derstanding the driving mechanism of mass-loss in AGB stars, 
one of the major open questions in stellar evolution. It is gen- 
erally accepted that the mass loss of AGB stars is dust-driven 
(Lamers & Cassinelli 1999). Pulsations bring the envelope far 
from the stellar equilibrium radius, where the temperatures are 
low enough for dust-condensation. This dust is easily acceler- 
ated by radiative pressure and it drags the gas along. For carbon- 
rich stars, the opacity of amorphous carbon is high enough in the 
near-infrared to drive mass loss. For oxygen-rich stars however, 
theoretical studies suggest that silicates are not able to drive the 
wind ( |Woitke|2006"l|Hofner & Andersen|2007l|Hofner|2008| l. It 
is clear that an exact understanding of the process is lacking. 

A recent study of the infrared spectra of nine S-type AGB 
stars from the Short Wavelength Spectrometer onboard of the 
Infrared Space Observatory (ISO/SWS), presented by |Hony| 
et al. ( 2009| l showed that the infrared spectra of S-type AGB 
stars can be significantly different from those of oxygen-rich 
AGB stars and already reported on the detection of dust emission 
at 28 fim attributed to magnesium sulfides. Although this study 
shows that S-type AGB stars exhibit several unique dust charac- 
teristics, the study is limited because ISO/SWS is only sensitive 
to the brightest (at infrared wavelengths) AGB stars and thus 
to the S-type stars with higher mass-loss rates, while a large and 
homogeneous sample is necessary to draw conclusions about the 
dust condensation in S-type stars. 

In this paper we present a combination of optical spec- 
troscopy, infrared photometry and infrared spectroscopy for a 
sample of 87 zVtfrom'cj^] S-type stars. This sample allows us to 
study the interaction between the chemistry and dynamics of the 
stellar atmosphere and the dusty circumstellar environment. In 
^2]and ^3]we discuss the sample selection procedure and the data 



reduction. In ^4] we derive the chemical and variability charac- 
teristics of the stars in our sample. In ^5] we give an overview 
of the "naked" stellar atmospheres and we show that their spec- 
tral appearance of depends strongly on the C/O ratio of the star. 
In ^6] we present a tool to decompose infrared (IR) spectra into 
the contribution of the stellar atmosphere and the different dust 
species, allowing us to classify the IR spectra of dust-rich stars 
into three distinct groups. We show that there is a relation be- 
tween the C/O ratio and those groups. 

2. The sample of S-type AGB stars 

In order to study the interaction between the stellar atmosphere 
and the infrared properties of S-type AGB stars, we selected a 
large sample of S-type stars designed to include all intrinsic S- 
type AGB stars bright enough to be detectable with the Spitzer 
IRS spectrograph. This sample spans the entire range of C/O 
ratios, temperatures and s-process abundances found in S-type 
AGB stars, but is slightly biased towards stars with lower mass- 
loss rates (see Sect. 1231). 



2 Intrinsic S-type stars are genuine thermally pulsating AGB stars, 
contrary to extrinsic S-type stars which are enriched in C or s-process 
elements by accretion from an evolved companion l Groenewegen |1993| 
Van Eck & Jorissenl 1999b. 



2.1. Sample selection 

A first sample was drawn from the second edition of the General 
Catalog of S stars (GCS S2, [Stephenson| 1 994) . This catalog clas- 
sifies 1346 stars as S-type stars, based on the presence of ZrO or 
LaO bands in the red part of the optical spectrum. This method 
only selects stars based on an enrichment of s-process elements 
in their lower envelope, which not only includes stars that have 
undergone a recent third dredge-up event, the intrinsic S stars, 
but also stars that are enhanced with s-process elements by bi- 
nary interaction, called extrinsic S stars ( Brown et al.|1990 i. In 
order to study the dust-sequence on S-type AGB stars, we need 
to exclude the extrinsic stars from our sample. 

Several properties distinguish both groups, such as lumi- 
nosity, the K - [12] color, the presence of technetium in the 
atmosphere, and variability ( |Van Eck & Jorissen| [l999). The 
best distinction method is based on the presence of technetium 
(Tc), which is brought to the surface during the third dredge-up 
events. With a laboratory half-life time of approximately 2 x 10 5 
years, the detection of Tc in the spectrum points to a recent third 
dredge-up event. Although the detection of Tc is the most se- 
cure way to include only intrinsic stars in the sample, it requires 
high resolution spectroscopy in the blue region where these cool 
stars are faint. A criterion based on the luminosity would in- 
troduce strong and unwanted biases on the distance or intrinsic 
brightness of the stars, while putting a lower limit on the infrared 
color K — [12] would favor stars with a dusty circumstellar en- 
vironment. Considering all this, source variability turns out to 
be the most efficient way to discriminate between intrinsic and 
extrinsic S-type stars. 

The AS AS light curves of the sample presented by |Van Eck| 
|et al.| ( [2000| were examined in the light of the intrinsic/extrinsic 
classification of these authors. For the extrinsic stars (excluding 
the symbiotic binaries), we found an upper limit of 0.12 mag 
for the standard deviations of the V-band light curves. Hence 
applying a lower limit of 0.12 mag on the standard deviations 
of the V-band light curves should result in the exclusion of the 
extrinsic S-type AGB stars. The remaining sample consisted of 
459 presumably intrinsic S-type stars with variability amplitudes 
larger than 0.12 mag. 

To estimate the source variability of the S-type stars, we 
have cross-identified all stars in the GCSS2 with the two ma- 
jor available variability surveys: the All Sky Automated Survey 
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(ASAS) loc ated at the Las Campanas Observatory (Pojmanski 
et al.||2005| l and the Northern Sky Variability S urve>p1 (NSVS) 



operated at Los Alamos National Laboratory ( Wozniak et al. 
2004a). This resulted in 1154 S-type stars with V-band light 
curves. 

In the last step, the sample was reduced to contain only stars 
that could be observed by Spitzer within a reasonable integration 
time. Therefore we used a cross-identification of the remain- 
ing sources with the catalogues of the Infrared Astronomical 
Satellite (IRAS) and the Two Micron All Sky Survey (2MASS). 
From the relation between the K magnitude and the 12 fim IRAS 
flujQwe estimate that S stars with 4.5 mag < K < 5.5 mag ex- 
hibit maximum fluxes between 0.3 and 3 Jy in the Spitzer wave- 
length range, not so bright as to saturate the detector but suffi- 
ciently bright to allow short integration times. The final target 
list contained 90 sources (3 of which will be rejected based on 
the ground-based optical spectra). 



2.2. Extrinsic intruders in our sample 

Because the cut-off criterion on the standard deviations of the 
V-band light curve is not strict, we can expect a small number of 
extrinsic S-type stars in our sample. As said in the previous sec- 
tion, the best way to (a posteriori) confirm that a star is an S-type 
AGB star is to check for the presence of the 42 38 A , 4262 A and 
4297 A Tc lines in high-resolution spectra ( Uttenthaler et al. 
[2007). 

As a filler program, we obtained high-resolution spectra for 
some of the brightest stars in our sample. The target selection 
was inhomogeneous, focussing on (i) stars with low-amplitude 
light curves and (ii) stars with an unusual appearance in the 
low-resolution optical spectra. The first group of stars was tar- 
geted because the chance of finding an extrinsic S-type star 
amongst the low-amplitude pulsators is larger since true AGB 
stars are expected to have larger pulsation amplitudes. The stars 
with an unusual low-resolution optical spectrum could be mis- 
classified as S-type stars in the GCSS2 (Stephenson 1994) (e.g. 
BD+02 4571). For more information on the low- and high- 
resolution spectra and the lightcurves, see Sect. [3] 

As shown by |Uttenthaler et aTT ( 2007[ l, we can make the dis- 
tinction between Tc-rich and Tc-poor stars by comparing the 
flux ratios in the Tc line to a pseudo-continuum flux. Using this 
method a reliable distinction of this kind can be made, even for 
low SNR spectra. 

Three stars did not show significant Tc lines and were re- 
moved from our sample (CPD-19 1672, BD+02 4571 and 
CSS 718). Because of the inhomogeneous subset of stars with 
high-resolution spectra, this ratio of 3 stars without Tc lines in 
a subset of 13 stars is not representative for the entire sample. 
Discarding these stars, we are left with a total sample of 87 S- 
type AGB stars. 



2.3. Bias towards low mass-loss rates 

During the target selection phase, we put a limit on the /f-band 
magnitude and cross-identified the sample with the ASAS and 
NSVS catalogs. Because both surveys are limited to V-band 
magnitudes of approximately 14 and the Zf-band magnitude is 



limited to the 4.5-5.5 mag range, we have effectively limited the 
sample to stars with V - K < 10 mag. 

Based on the marcs model atmospheres for S-type AGB stars 
without detectable circumstellar material (see Section|4|, we can 
expect V - K values up to 10 mag or more for stars with 7^ 
lower than 2700 K, even without circumstellar material. Stars 
with circumstellar dust are reddened and thus have higher val- 
ues for the V - K color. The limitation on this color will exclude 
the stars with the highest mass-loss rates because they appear 
redder ( jGuandalini|2010| >. 

In Fig. [1] we show the J — K and V - K colors for our stars, 
the sample of S stars observed with ISO/S WS ( |Hony et al.|2009") 
and all S-type stars found in the Stephenson Catalog for compar- 
ison. For our stars, the phase-calibrated V-band magnitudes are 
derived from the light curves in the ASAS or AAVSO database 
and new measurements with the Mkll photometer, discussed in 
more detail in Sect. |3.2.1| For the ISO/S WS sources we use the 
AAVSO database to derive the phase-calibrated V magnitudes. 



3 See http://skydot.lanl.gov 

4 Actually, the IRAS fluxes are given in Jy and are thus monochro- 
matic flux densities, but for ease of use, we will continue to use the term 
flux. 



The relations presented in Carpenter (2001 1 were used to con- 
vert the Admass magnitudes to -Ksaao magnitudes. For the S- 
type stars found in the Stephenson catalog that are not observed 
with Spitzer or ISO/SWS, the V - K and J - K colors are not 
phase-calibrated and not corrected for the interstellar reddening 
because the distances and pulsation characteristics are very un- 
certain. 

From this figure, we can see that although we have indi- 
rectly removed stars with V — K > 10 mag from our sample, 
this effect is only small and we excluded only the coolest stars 
which have high mass-loss rates, less than 5% of all S-type stars. 
Furthermore, the figure shows that the sample of S-type AGB 
presented in Hony et al. ( 2009) > is biased towards higher mass- 
loss rates, because the hmited sensitivity of the ISO/SWS spec- 
trograph favors stars with a strong infrared excess. We have to 
keep both biases in mind when we compare our results with lit- 
erature. 



3. Data 

3.1. Spitzer Space Telescope 

The IRS data of all 90 S-type stars were obtained using the IRS 
Standard Staring mode. The targets were observed in two posi- 
tions in each aperture (SL1, SL2, LL1 and LL2) to obtain low- 
resolution spectra over the entire wavelength range (5-38 //m). 
All data were taken during a period from September 2006 to 
October 2007 (Program ID 30737, P.I. S. Hony). 

For the data reduction we used the FEPS pipeline, devel- 
oped for the Spitzer legacy program "Formation and evolution 
of planetary systems". A detailed description can be found in 
|Hines et al.| ( |2005] ). For the extraction we start from the inter- 
mediate droop data products as delivered by the Spitzer Science 
Center, together with the SMART reduction package tools de- 
scribed in detail in Higdon et al. (2004). The spectra were ex- 
tracted using a fixed width aperture, so that 99% of the source 
flux falls within the aperture. After the extraction of the spec- 
trum for each nod position and cycle, a mean spectrum over all 
slit positions and cycles is computed for each individual order. 
Thereafter, the orders are combined. In regions where there is 
order overlap, the fluxes are replaced by the mean flux value 
at each wavelength point. The IRSFRINGE package was used 
for the de-fringing of the spectra. For the calibration, spectral 
response functions were calculated from standard stars and cor- 
responding stellar models. The estimated average uncertainty on 
the spectral response function and uncertainty on the extracted 
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Fig. 1. The V — K vs. J - K diagram of the S-type stars in the 
Stephenson catalog (light grey dots), the ISO/SWS targets pre- 
sented in Hony et al. |2009) (black diamonds) and the Spitzer 
sample presented in this paper (dark grey squares). 



spectra are incorporated in the final estimated error. Finally, for 
each target a one-dimensional spectrum is computed by stitching 
together overlapping orders. This is done by multiplying the LL 
orders with a correction factor to match the mean flux with the 
SL orders in the overlapping wavelength range. No wavelength 
shifting or spectral tilting is performed. 

3.2. Additional ground-based data 

3.2.1. SAAO 

For 69 of the 87 stars in our sample, we have obtained 
JHKL magnitudes (1.25, 1.65, 2.2 and 3.5 //m) at the South 
African Astronomical Observatory (SAAO). The observations 
were made with the Mkll infrared photometer on the 0.75-m 
telescope. The accuracy is typically better than 0.03 mag in the 
J, H, and K band, and better than 0.05 mag in the L band. All 
JHKL magnitudes derived for these stars are given in Table [T] 

3.2.2. Low-resolution optical spectra 

Low resolution optical spectra were obtained for 81 of the 87 
stars in our sample. The southern objects were observed with 
the 1.9-m telescope at SAAO, whereas the more northern ob- 
jects were observed with the ISIS spectrograph mounted on the 
4.2-m William Herschel Telescope (WHT) at the Observatory de 
Roque de los Muchachos (ORM) situated at La Palma. 

The spectra were wavelength calibrated with a CuNe+CuAr 
lamp. The correction for atmospheric extinction was performed 
using the average curve for the local, continuous atmospheric ex- 
tinction at the site. The resulting SAAO spectra cover the 4800- 
7700 A range with a resolving power of R = A ~ 3800. The 
spectra of WHT consist of a blue arm, spanning 3700-5350 A 



with a resolving power of R » 5000 and a red arm, spanning 
5100-8100 A with a resolving power of R » 3800. 

3.2.3. High-resolution optical spectra 

We obtained high-resolution spectra for some of the brightest 
targets in our sample. Eight southern sources were observed 
with the Coralie echelle spectrograph at the Euler telescope 
in La Silla, spaning 3900-6800 A, with a resolving power of 
R ~ 50 000. Thirteen northern sources were observed with the 
HERM ES echelle spectr ograph at the Mercator telescope in La 



Palma (Raskin et al. 201 1 1, covering 3770-9000 A with a resolv- 
ing power of R « 85 000. We use the end-products of the auto- 
matic data reduction pipelines for both instruments. 



3.2.4. Dereddening 

To remove the effects of interstellar, wavelength dependent ex- 
tinction from the ground-based and Spitzer data, we start from 
the 3-dimensional reddening models of Drimmel et al. (2003) 
and |Marshall et al.| ( |2006) . The Marshall model returns the K- 
band extinction for over 64 000 lines of sight, each separated 
by 15', based on the 2MASS point source catalogue. When 
possible, we use this model to estimate the extinction, but the 
Marshall model only covers the inner Galaxy (|/| < 100°, \b\ < 
10°). The model presented in |Drimmel et al. (20031 returns V- 
band extinctions for the entire sky. The data are dereddened, us- 
ing the extinction laws of Cardel lTet al.| ( |1989| l for the optical 
and NIR data and |Chiar & Tie lens ( 2006jTat longer wavelengths. 



4. Stellar parameters 

Because the spectral appearance of S-type AGB stars depends 
strongly on the stellar parameters such as T e ff, C/O ratio, s- 
process element abundance and iron abundance, it is possi- 
ble to constrain these stellar parameters by comparing model 
atmospheres to observational data. These stellar parameters 
characterize the chemical conditions in the stellar atmosphere. 
Especially the C/O ratio is expected to have a very profound im- 
pact. 

Because the sample selection is based on the ASAS and 
NSVS catalogs, we are able to constrain the pulsation charac- 
teristics of most stars. The period and amplitude of the V-band 
variations characterize the variability of the stellar atmosphere. 

4.1. Chemical characterization 

4.1.1. Comparison with MARCS model atmospheres 

The stellar parameters were derived using a grid of model atmo- 
spheres specially computed for S-type stars with the marcs code 
( Gustafsson et al.|2008| l. This model grid has been described in 



|VanEck et al (2010 



and Neyskens et al. (2010!. The full de 



tailed description will be published in a forthcoming paper (Plez 
et al., in preperation). To summarize briefly, the grid covers the 
following parameter space: 

- 2700 < r eff (K) < 4000 in steps of 100 K 

- C/O = 0.5, 0.750, 0.899, 0.925, 0.951, 0.971, 0.991 

- [s/Fe] = 0, +1, +2dex 

- [Fe/H] = -0.5, dex. 

The unequal C/O spacing has been designed to cover in an opti- 
mal way the change of the absorption features affecting the op- 
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tical spectra when C/O increases from 0.5 to unity. All mod- 
els were computed for a stellar mass of M = 1 M and with 
[a/Fe] = -0.4 x [Fe/H]. It is essential to take into account the 
non-solar C/O ratio of S-type stars since it has an important im- 
pact on the partial pressures of H2O and TiO, two major opacity 
contributors. Coupled with s-process overabundances, it affects 
the depth of the prominent ZrO bands observed in S-type stars. 
The effect of the gravity was found to be too degenerate to be 
constrained solely with the observational data available for the 
present study. Therefore it was fixed in the following way, de- 
vised to approximately follow the portion of the AGB covered 
by the S stars under consideration: log g = if T e g < 3400 K, 
logg = 1 otherwise. Similarly, low-resolution spectra do not 
allow to lift the degeneracy between metallicity and s-process 
overabundance. 

A set of well-chosen photometric and narrow-band indices 
were then computed on the basis of low-resolution synthetic and 
observed spectra (see Van Ecket al. ( 2010 1 for more details). The 
low-resolution spectra used here consist of Boiler & Chivens 
spectra (AA = 3 A, 4400 - 8200 A), SAAO spectra (A^ = 1.8 A, 
4800 - 7700 A) and WHT spectra (AA = 1.8 A, 5100 - 8100 A). 
Two colour indices V-K and J-K as well as three spectroscopic 
indices (strengths of the carefully-selected ZrO and TiO bands, 
as well as the strength of the Na D lines) were used to define a 
total of five indices. 

Chi-square minimization between observed and synthetic in- 
dices then led to a set of plausible physical parameters (effec- 
tive temperature, gravity, metallicity, C/O ratio, and [s/Fe]). The 
agreement between selected synthetic spectra and the observed 
optical and infrared low-resolution spectra is in most cases very 
good. However, it is sometimes difficult to select one out of the 
first few "best-matching" synthetic spectra. In fact, the uncer- 
tainties inherent to the whole procedure (e.g., non-simultaneous 
photometry and spectroscopy, stellar variability, or approximate 
values of the reddening) have to be taken into account in estimat- 
ing the best model and the uncertainty on the stellar parameters. 
The resulting value and uncertainties on the stellar parameters 
r e ff , C/O, and [s/Fe] that provide a reasonably good fit between 
the low-resolution observed and synthetic spectra are given in 
Table [3] 



4.1.2. Stellar properties 

From the results in Table [3] we can get a general overview of the 
chemical properties of the stars in our sample. The effective tem- 
perature spans a range from 2900 K for the coolest S-type star 
to 3650 K for the hottest one, with most stars with T e ff around 
3300 K. There is a strong correlation between the place in the 
ZrO-TiO plot and the C/O ratio. The ZrO and TiO indices can 
thus be used as indirect indicators for the C/O ratio (see Fig. [2J. 
We will make use of this relation throughout the entire paper. 
The figure shows that one can separate two groups, stars with a 
low C/O ratio are located at the upper left corner and the stars 
with a high C/O ratio towards the lower right corner. The posi- 
tion in this diagram can thus be used as a decent estimate for the 
C/O ratio. 

From this figure, we can also see that there is a significant 
correlation between the [s/Fe] values and the C/O ratios. This ob- 
vious correlation was predicted and observationally confirmed in 
AGB stars ( [Smith & Lambert|1990||Mowlavi|1997l l and S-type 
stars in particular ( Van Eck et al.||2010[ ). Because all combina- 
tions of C/O and [s/Fe] were allowed, the presence of the relation 
is not inherent to the method to derive the stellar parameters, but 
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Fig. 2. The average line-to-continuum-ratios of the ZrO and TiO 
absorption bands within the spectral range of WHT and SAAO 
low-resolution spectra (i.e. ZrO and TiO indices) for all stars in 
the sample. The size of the symbols increases with [s/Fe] and the 
color varies with C/O. Crosses indicate stars for which T ei f, C/O 
and [s/Fe] are unknown. 



rather shows that the method works well and can reproduce this 
relationship. 



4.2. Variability characterization 

To derive the pulsation parameters, we use the databases of 
NSVS and AS AS. A complete description of the instruments 
and reduction processes can be found in Wozniak et al. ( 2004b) 
and [Pojmariski (2004). For some stars, additional V-band light 
curves are available from the American Association of Variable 
Star Observers (AAVSO). 

Because for most stars only a few data-points are available, 
we use the generalized least-squares periodogram to find the 
pulsation periods of all stars (Zechmeister & Kurster 2009 ). This 
periodogram takes into account that for small datasets the sam- 
ple average can differ significantly from the true average. Some 
semi-regular variables are known to exhibit long secondary peri- 
ods (LSP), typically 5 to 15 times longer than the primary period. 
When detected, such variations are removed (prewhitened) from 
the light curve and the period analysis is repeated. The results of 
this analysis can be found in Table [2] 

The period analysis resulted in detected periods for 54 of the 
87 stars in our sample. For most cases where no significant pe- 
riodicity was found, this can be explained by the short timespan 
of the NSVS light curves. These short timespans allow an iden- 
tification of the amplitude of the variations, but a period analysis 
needs longer timespans. For eight stars, the lack of periodicity is 
due to intrinsic irregular behavior. 

In addition to these periods, the General Catalog of Variable 
Stars (GCVS) gives the period to 18 stars. For 11 of these 
stars, we also derived a period, based on the NSVS, ASAS and 
AAVSO light curves. As shown in Table|2j the derived period for 
10 of these 1 1 stars is in close agreement with the periods shown 
in the GCVS. However, the recent data shows that V899 Aql has 
a period of 374 days and is most likely a Mira, this star is clas- 
sified as a semi-regular variable in the GCVS, with a period of 
100 days. In further calculations, we will use the periods and pul- 
sation types derived from the NSVS, ASAS, and AAVSO light 
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curves if available. The sample contains 6 irregular pulsators, 38 
semi-regular variables and 21 Miras. 



5. Naked stellar atmospheres 

The 87 S-type stars in our sample show a wide variety of features 
in the infrared spectra. Some stars show clear dust/molecular 
emission features, while others only show the molecular absorp- 
tion features typical for "naked" stellar atmospheres. Because 
stars with and without emission features have to be treated sepa- 
rately, we classify the targets into naked stellar atmospheres and 
stars with additional emission features. Out of the 87 stars, 47 
are classified as naked stellar atmospheres. The other stars show 
relatively strong emission features on top of these stellar atmo- 
spheres which are due to (i) molecular emission from extended 
atmospheres (3 stars), (ii) the typical emission features from hy- 
drocarbons (4 stars) or (iii) dust emission (32 stars). Here, we 
focus on the naked stellar atmospheres. The dust emission will 
be discussed in more detail in Sect.[6]and the molecular emission 
in Sect. 13 



5. 1 . Overview of the observed absorption bands 

In the bottom panel of Fig. [3] we show the normalized infrared 
spectra of all 47 naked stellar atmospheres of our sample; the 
top panel shows the normalized templates used for compari- 
son. These templates have column densities and temperatures 
that are representative for AGB stars (CO at T = 2000 K and 
N = 10 21 cm 2 ; H z O at T = 1500 K and N = 10 20 cirT 2 ; SiO 
at T = 2000 K and N = 10 21 cm" 2 ; SiS at T = 1500 K and 
Cami|2002l|Cami et al.|2009| >. 



N 



10 19 cm" 2 ; 



From this figure it is clear that all stars in our sample show 
absorption bands of oxygen-rich molecules such as H2O and 
SiO. The right panel shows the 19 stars with a significant SiS ab- 
sorption band at 13/zm. The remaining panels show the 28 stars 
without SiS absorption. We did not detect the strong absorption 
features at 7.2 (S0 2 ), 7.6 (H 2 2 ), 13.8 (C 2 H 2 ), 14 (HCN) or 
15 //m (CO2) typical for of S0 2 , H 2 2 , 2h2, HCN and C0 2 . 



5.2. The relation between SiS absorption and the C/O ratio 



Considering the partial pressures, |Cami et al. (2009) predicted 
that SiS absorption will only be visible in the S-type AGB stars 
with C/O ratios very close to unity. Fig.|4]shows the ZrO and TiO 
index for all naked stars, with a color scale indicating the C/O 
ratio for stars with and without significant SiS absorption. From 
this figure it is clear that there is a strong relation between the 
estimated C/O ratio and the presence of SiS absorption: all stars 
with C/O > 0.96 show SiS absorption and all stars with C/O < 
0.96 have no significant SiS absorption bands. It should be noted 
that the C/O and SiS determinations are completely independent. 
The C/O ratio is determined, based on the optical spectroscopy 
and optical and NIR photometry, while the SiS bands are only 
observed in the Spitzer spectra. 

Two conclusions can be drawn from this: (i) the C/O ratio 
is key to understanding and modeling the IR spectra of S-type 
AGB stars and (ii) the presence of the SiS absorption bands can 
help constrain the C/O ratio in stars, which can be a useful tool 
when the C/O ratio cannot be derived otherwise. 
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Fig. 4. The ZrO and TiO indices for the naked stellar atmo- 
spheres in the sample. Squares and circles indicate stars with 
and without significant SiS absorption with a colorscale that in- 
dicates the C/O ratio. Stars for which we have no estimates for 
T e ff, C/O, and [s/Fe] are indicated with open symbols. The figure 
shows a one-one relation between C/O ratio and presence of SiS 
absorption in the infrared spectrum. 



6. Dust identification 

In this section we focus on the 32 stars with significant dust 
emission features in the IR spectra. We first give an overview 
of the dust formation sequence in M- and C-type stars and then 
we compare this with what we find for S-type AGB stars. 



6.1. Dust formation in AGB stars 

The study of the dust condensation chemistry in AGB stars has 
mostly focused on the formation of oxygen-rich dust species in 
M-type stars or carbon-rich dust species in C-stars. In Fig.|5]we 
show the formation pathways for several dust species that are 
known to be important in oxygen-rich and/or carbon-rich stellar 
outflows. 

Although it is not yet observationally confirmed, a consis- 
tent theory for the dust formation in M-type AGB stars exists. 
Theoretically, the first condensate to form is alumina grains 
(Al 2 03), which can condense at temperatures as high as 1760 K. 
When these grains interact with the circumstellar SiO and Ca, 
melilite is formed. This melilite is first pure gehlenite, but can 
later on become akermanite due to additional interactions with 
the circumstellar Mg. Due to the formation of akermanite from 
gehlenite grains, alumina grains are released from the gehlenite 
grains and can go on to form additional alumina grains (if the lo- 
cal temperature is above 1510K) or spinel grains (if the temper- 
ature is below 1510K). The melilite reacts at 1450 K to produce 
diopside and more spinel. The reaction of diopside and spinel 
can form anorthite. Further on in the wind, at lower tempera- 
tures, these grains can be the seeds for the formation of enstatite, 
forsterite and the amorphous silicates we observe in many AGB 
stars, (this scenario was proposed and adapted by Grossman & 
|Larimer|[T974l |Onaka et al.|[T989l |Tielens|[l990l |cami||2002) . 
Additionally, several species, that are not part of the formation 
path that eventually ends with silicates, can condense at certain 
temperatures. For example, due to the reaction of H 2 with Fe, 
iron oxides form at a temperatue of approximately 900 K, while 
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Fig. 3. Top panel: the template spectra of SiS, SiO, H2O, and CO. Bottom panel: the spectra of all naked stars from 5 pm to 20 //m, 
normalized with a blackbody spectrum. The the two panels on the left show the stars without SiS, the right-most panel shows all 
stars with significant SiS absorption. The sharp feature at 15.3 /im seen in all stars is most likely an artifact of the data reduction 
process. 



the magnesium-rich oxides can form at higher temperatures (up 
to 1200 K) due to the reaction of H2O with atomic magnesium 
( |Ferrarotti & Gail|2003| >. 

The broad picture for C-type AGB stars is simpler, but some 
key details have yet to be filled in. The main dust species de- 
tected in carbon-rich stars are amorphous carbon, silicon car- 
bide and magne sium sulfide (jTreffers & Cohen|1974 Goebel & 



Moseley 1985 



Chan & Kwok 



Martin & Roge rs||1987 Andersen et al. 
1990| JHony et al.||2002[ |Zijlstra et al. 



1999 



2006 



Lagadec et al. 2007 Leisenring et al. 20b8[ ). One of the key 
questions in the carbon-species formation is the importance of 
mantle chemistry, for example the nucleation of amorphous car- 
bon and magnesium sulfides as mantles around silicon carbide 
grains have been proposed as important dust formation mecha- 
nisms (Lorenz -Martins et al. 2001 1 Zhukovska & Gail 2008 ). 

Assuming a certain dust condensation sequence and an ex- 
panding stellar wind, one expects a process called freeze-out, 
where the sequence is not completed but stops at an intermediate 
product when the wind density drops below the density required 
for the next step. This freeze-out has been tentatively observed 
in AGB stars (e.g. |McCabe et al.|1979||Heras & Hony|2005[ ). An 
alternative, but qualitatively similar effect could be the exhaus- 



tion of oxygen in these S-type environments. If the C/O ratio is 
close to unity, the formation of alumina grains can use up most 
of the oxygen atoms available after the formation of CO. If this 
happens, the density of free oxygen atoms is too low to continue 
the dust condensation sequence along the path presented in Fig. 
[5] ana logous to the freeze-out (as predicted by Sloan & Price 
1998). We will refer to this process as oxygen depletion. 



6.2. Dust decomposition tool 

The first crucial step towards understanding the dust condensa- 
tion sequence of these S-type stars in the present sample is the 
identification of the dust emission features in the infrared spec- 
tra. We developed a simple tool to decompose the infrared spec- 
tra into a stellar continuum and the dust emission features of 
each individual dust species. 

If we assume an optically thin circumstellar environment, 
where the dust condenses out instantaneously once the tempera- 
ture gets below a certain threshold Tq and where the dust grains 
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Fig. 5. The estimated dust condensation temperatures for different dust species found in AGB stars. The large inset panel zooms 
in on the dust condensation sequence for oxygen-rich stars, starting from alumina grains. We highlighted the species which were 
included in our dust decomposition tool in grey. 

are in thermal equilibrium, we can write the temperature profile where rn corresponds to the radius where the temperature in the 
in the circumstellar layers as: outflow drops below Tq. 



TdustO) = r< 



■1/2 



The flux coming from a single dust species (using spherical 
' grains and Mie theory) in an optically thin stellar wind with a 



8 



K. Smolders et al.: The Spitzer Spectroscopic Survey of S-type Stars 



constant mass-loss rate M4 and outflow velocity v can be written 

as: 



v, dust 



3M d & ro_ 
4vp d a D 2 



B v (T d (x))dx 



(2) 



where is the density, Q Y the absorption efficiency, a the radius 
of the spherical dust grains, ro the inner radius of the dust shell, 
D the distance to the object and T&(x) the temperature profile as 
in Eq.[T]at a distance x - r/ro ( Schutte & Tielens 1989 1. Finally, 
we can combine all constants and unknown parameters into a 
single parameter B, for each dust species ;. The total spectrum 
can then be written as: 



Q vJ B v {T A (x))Ax 



(3) 



where F Vy * is the flux of the central star at frequency v and B, 
is the scaling constant that controls the relative strength of the 
emission of a specific type of dust. 

To arrive at the best fitting model, we use the Lawson and 
Hanson Non-Negative Least-Squares algorithm (NNLS), which 
gives the best model under the condition that B, > 0. To es- 
timate the uncertainties on the B, parameters, we multiply the 
uncertainties on the data with a factor until the x 1 value for the 
best fit is 1 . Using these new uncertainties, we add normally dis- 
tributed noise to the spectrum and calculate the scaling factors 
Pi for 1000 randomly disturbed spectra. If the mean value is sig- 
nificantly non-zero (i.e. outside the 99.9% confidence interval), 
we consider this a positive detection of that dust species in the 
infrared spectrum. 

Because the MARCS model atmospheres still need valida- 
tion at infrared wavelengths, we prefer to fit the stellar contin- 
uum and the emission simultaneously. Therefore, we allow the 
dust decomposition tool to make any linear combination of the 
naked stars presented in Sect. [5] and Planck functions of differ- 
ent temperatures, ranging from 500 K to 4000 K, as an estimate 
for the stellar continuum, F Vt *. The big advantage of using the 
naked stars is that we can change the relative depth of the CO, 
SiO, H2O, and SiS bands to give the optimal fit, while we can 
change the slope of each individual star with the Planck func- 
tions. The drawback, however, is that we will not be able to dis- 
tinguish between the stellar continuum and additional feature- 
less continuum emission from e.g. metallic iron or amorphous 
carbon. 

Finally, we have to decide which dust species to include in 
the dust decomposition tool and which to leave out. This is of 
vital importance because too few dust species will leave us un- 
able to explain some emission features, while too many species 
would make the interpretation difficult and would lead to degen- 
eracy. We started from a wide range of dust-species, including 
titanium oxides, titanium carbides, silica and all dust species pre- 
sented in Fig. [5] We rejected those species that did not improve 
the goodness-of-fit to the IRS spectrum. The final selection of 
dust species is: 



Amorphous silicates The strongest emission features due to 
oxygen-rich dust are the two emission features at 9.8 and 
18yt/m due to the stretching and bending of Si-O bonds in 
amorphous silicates. These features dominate the infrared 
spectra of most M-type stars. Amorphous silicates form at 
relatively low temperatures and are thus expected to be at 
the end of the dust condensation sequence. Here, we have 
chosen olivines as representative of the typical silicates that 
are observed in AGB stars. Optical constants have been taken 
from |Jaeger et al.| ( [T994]). 

Gehlenite We know from Hony et al. ( 2009 \ that the peak of 
the silicate features in S stars is not at 9.8 fim as in M-type 
stars, but shifted towards redder wavelengths, up to 10.6 fim. 
Since this cannot be explained with only amorphous sili- 
cates, we need to include gehlenite, which is an alumino- 
silicate (Ca2Ai2Si07). A combination of amorphous sili- 
cates, gehlenite and amorphous alumina can explain the 
broad and shifted 10.6 micron feature. Optical constants 
have been taken from IMutschke et aTT] < fT998> . 

19.8 fan feature Although the identification of this distinct 
emission feature at 19.8 fim is often attributed to iron ox- 
ides (Cami 2002; Posc h et al.| 12002), other authors have 
argued that the 19.8 fim emission, together with additional 
emission features at 28 and 32 fim is due to crystalline sili- 
cates (Slo an et al.|2003 1. These features are visible in many 
M-type stars and most clearly in stars with low mass-loss 
rates. We were unable to reproduce the 19.8 fim emission us- 
ing only crystalline silicates without introducing strong addi- 
tional emission features in the 9-1 1 fim region, which we do 
not observe. Therefore, we prefer iron oxides (more specifi- 
cally Mgo.1Feo.9O). Optical constants have been taken from 



Henning et al. ( 1995 1. 



The 13 fim feature Since the first detection of the 13 fim fea- 



carriers: corundum (crysta 



ture by Vardya et al. ( 1986 1, it has been attributed to different 



lline AI2O3 , [Glaccum| 1995 , silic a 
(Si0 2 , |Speck|l998) and spinel (MgAl 2 a ^Posch et al.|1999) . 
The debate on the exact carrier of the 13 /mi feature is still 
ongoing (DePew et al. 2006 1. In our dust decomposition tool, 
we tried all three carriers and found that both corundum and 
spinel can reproduce the emission, but the best results are ob- 
tained with spinel. However, detailed modeling is required to 
attribute the 13 fim unambiguously to any one carrier, which 
is beyond the scope of this paper. Optical constants for spinel 
have been taken from Fabian et al. (2001). 
Magnesium sulfides The most remarkable feature in the spec- 
tra of these S-type stars is the broad 26-30/im emission fea- 
ture, attributed to magnesium sulfides. A similar feature has 
already been observed in two S-type stars, W Aql and n l Gru 
( |Hony et al. 2009). Because of the lower spectral resolu- 
tion of our set of spitzer spectra, the two separate emission 
peaks at 26 and 28 fim visible in the ISO/SWS spectra are 
blended to a single, broad emission feature. In our dust de- 
composition tool we use the mass absorption coefficients of 
magnesium sulfide. Optical constants have been taken from 



Begemann et al. ( 1994 ). 



Amorphous alumina Amorphous alumina has a high conden- 
sation temperature and is thus expected to be at the start of 
the dust condensation sequence of oxygen-rich stars (see Fig. 
[5}. The most prominent feature of amorphous alumina is at 
11.2/jm with a broad red wing, reaching up to 15.6 /im. It 
is a common component of the circumstellar wind of M- 
type AGB stars. Optical constants have been taken from 
IBegemann et al.| ( fT^7) . 



Although we could expect the presence of the narrow 
113 fim emission feature of amorphous SiC grains (optical 
constants have been taken from Mut schke et al. |1999| ) in stars 
with C/O ratios close to one, replacing one of the aforemen- 
tioned dust species with amorphous SiC grains resulted in higher 
X 2 values for all spectra. We can confidently state that amor- 
phous silicates, gehlenite and amorphous alumina are necessary 
to explain the data. Furthermore, the inclusion of SiC, together 
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with the other dust species, did not significantly lower the;^ 2 for 
our sample of S-type stars. Therefore we did not include amor- 
phous SiC grains in the simple fitting routine. More detailed 
modeling is necessary to make a statement on the absence or 
presence of SiC grains, but we can conclude that the SiC emis- 
sion is not strong in any of our stars. 

6.3. Overview of the detected dust species 

The results of the dust decomposition tool can be found in the 
online Appendix. Although S-type stars can have C/O ratios very 
close to unity, the majority of the stars only show dust emission 
features similar to those found in oxygen-rich stars. 

6.3.1 . Classification of the dust spectra 

Based on the results of the dust decomposition tool, we can 
subdivide the infrared spectra into 3 groups of spectra with: 
(i) strong emission features of silicates and gehlenite at 10 and 
18 /mi (16 stars), (ii) emission of alumina at 11. 2 /mi and of 
magnesium sulfide at 26-30 //m (10 stars), or (iii) a 13 /mi and 
19.8 //m emission feature (6 stars). 

Group I, stars with strong silicate/gehlenite features: 

Approximately half of the stars with dust features in our 
sample, show a dominant emission feature around 10/mi, 
due to the stretching of Si-O bonds of silicate or gehlenite 
dust grains. In Fig. [6] we show all stars with a strong 10 /mi 
feature. All spectra in the right panel show substructure at 
9.8 and 1 1 /mi, sometimes accompanied by a sharp peak 
at 13 /mi, the spectra in the left panel show a single broad 
emission feature that peaks at wavelengths up to 10.9 //m. 
Although the dust decomposition tool can reproduce the 
emission features with substructure, it is impossible to 
reproduce the broad and smooth emission features, peaking 
at 10.5 //m, visible in the left panel, without showing 
substructure. 

|Hony et al. ( 2009 ) found 1 1 stars with a single smooth emis- 
sion feature at ~ 10.5 //m. This feature was only found in S- 
type AGB stars, while there were no oxygen-rich AGB stars 
that exhibit the same emission feature. The authors argue that 
the broad and smooth emission found in these S-type stars is 
fundamentally different from the substructured emission in 
oxygen-rich AGB stars and the shift towards longer wave- 
lengths is due to an increase in the magnesium-to-iron ratio 
in S-type stars. 

Group II, stars with alumina, 13 /<m and 20 //m features: 

Fig. [JJ shows the six stars in Group II. Although some stars 
do show some indication for the presence of silicates, the 
dominant emission features are the 11.2 /mi emission due 
to amorphous alumina, a strong and sharp 13 /mi feature 
and the broad 19.8/mi emission. A similar combination is 
also common in oxygen-rich stars with low mass-loss rates 
dSloan et al.|2003] >. 

Many stars with a strong 13 /mi feature also show the pre- 
viously detected, unidentified emission features at 28 and 
32 //m indicated with arrows in Fig. [JJ Because the 28 and 
32 /mi features can blend into a broad plateau with two nar- 
row emission features at 28 and 32 /mi, the dust decomposi- 
tion tool sometimes wrongly models this broad plateau with 
emission of magnesium sulfides. This error can easily be 
identified because the residuals still show the narrow 28 and 
32 /mi emission peaks. 
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Fig. 6. The spectra of all stars in our sample with an infrared 
spectrum dominated by the 10 and 18 /mi features of silicates 
and gehlenite (Group I). The top left panel shows the smoothest 
features at 10 /mi, while the features in the top right panel show 
substructure at 9.8, 1 1.2 and/or 13 /mi. The bottom panel shows 
CD-39 2449, with the relative flux contribution of silicates and 
gehlenite. 



Several carriers have been suggested for these features, 
the most important candidates being crystalline silicates. 
However, the dust decomposition tool cannot explain the 28 
and 32 /mi emission features with these carriers without in- 
troducing strong additional emission features that we do not 
observe. 

Group III, stars with magnesium sulfide and alumina 
features: The ten stars in Fig. [8] show two dominant emission 
features, one attributed to amorphous alumina and the other 
to magnesium sulfide. Because of the absence of strong sili- 
cate features or a strong 19.8 /mi feature, the 26-30 //m emis- 
sion of magnesium sulfides is clearly visible on top of the 
11. 2 /mi of amorphous alumina and can easily be separated 
from the 28 and 32 //m emission of Group II stars. 
The substructure that can be seen in some of these stars at 
6.6 and 13.6 /mi is most likely due to emission of SiS gas 
near the star, as explained in Sect. 7.2 
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Fig. 7. The spectra of all stars in our sample with an infrared 
spectrum dominated by the emission of amorphous alumina at 
11.2/zm, the 13/mi and the 19.8/mi feature (Group II). The 
bottom panel shows CSS 480, with the relative flux contribu- 
tion each of these emission features. The previously detected, 
unidentified emission features at 28 and 32 //m are indicated with 
black arrows. 



6.3.2. C/O dichotomy 

In Fig.[9]we show the ZrO and TiO indices of all stars in our sam- 
ple with dust emission features, using different symbols for the 
different dust groups. From this figure it is clear that (i) the stars 
in Group I are distributed more or less uniformly throughout this 
graph, spanning a wide range of ZrO indices, TiO indices and 
C/O ratios, while (ii) the stars from Group II are located in the 
upper left corner of the diagram with low C/O ratios, and (iii) the 
stars from Group III clutter in the lower right corner with C/O 
ratios close to unity. 

Because the number of stars with known C/O ratios in each 
separate group is small, we used an unpaired T-test to check 
whether there is a significant difference in C/O ratios between 
these groups. We found that the C/O ratios of stars in Group 
II are significantly lower than the values for Group III stars 
(P < 0.1%). This dichotomy shows that there is a clear depen- 
dence of the circumstellar dust mineralogy on the chemical com- 
position of the stellar atmosphere in general and more specifi- 
cally on the C/O ratio. 



Fig. 8. The spectra of all stars in our sample with an infrared 
spectrum dominated by the emission of amorphous alumina at 
11. 2 /an and of magnesium sulfides at 26-30//m (Group III). 
The bottom panel shows CSS 987, with the relative flux con- 
tribution of alumina and magnesium sulfides. 



6.3.3. Comparison with the Sloan & Price classification 



Sloan & Price ( 1995] 1998| l also present a classification scheme 



for IR spectra of AGB stars, using flux ratios at 10, 11, and 
12 /mi to quantify the shape of the spectrum in that region. The 
ratio Fio/Fi i is an a proxy for the width of the 10 yum feature and 
F10/F12 for the strength of the red shoulder. |Sloan & Price| (l995 ) 
find that these flux ratios are strongly correlated and represent 
this correlation with a powerlaw. Along this sequence, they sub- 
divide the stars into eight separate classes, where the stars with 
the lowest F^/Fn and F^/F^ ratios fall in the SE1 class and the 
stars with highest ratios in the SE8 class. 

The interpretation of this sequence is that the spectral ap- 
pearance of the stars in classes SE1-SE3 is dominated by oxides, 
the stars in classes SE6-SE8 are dominated by amorphous sili- 
cates and the intermediate classes show indications of both dust 
species. This change in spectral appearance with increasing SE- 
class number has been explained as either an effect of increas- 



ing mass-loss rates (for example in |Sloan et aL]|2003[ |Pitman 
|et al.|2010] > or an effect of decreasing C/O ratio ( |Sloan & Price 
1998). Both interpretations are based on the fact that alumina 
grains condense before silicate grains. If the binding of oxygen 
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Fig. 9. The ZrO and TiO indices for the stars with different types 
of dust emission features. Group I, II and III are respectively 
indicated with dark grey diamonds, black triangles and light grey 
squares. The stars without dust emission features are shown as 
dots for comparison. 
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Fig. 10. The classification scheme for IR spectra of AGB stars, 
presented in |Sloan & Rnc^ ( [T995l[T998| ). Group I, II and III are 
respectively indicated with grey diamonds, black triangles and 
light grey squares. The solid line shows the powerlaw presented 
in the original paper. The dashed lines indicate the boundaries 
between the subsequent classes of the Sloan & Price classifica- 
tion. 



and alumina uses up the available oxygen, no free atoms will re- 
main to form silicate grains, and alumina grains will dominate 



the resulting shell ( Stencel et al.|19 90). The difference between 
both interpretations is only whether the mass-loss rate or the C/O 
ratio is the dominating cause of this observed sequence. 

In Fig. 10 we compare this power law with the flux ratios 
we find for the stars in our sample of S-type stars. Although we 
subtract the continuum flux in a different way, using the contin- 
uum estimates determined by the dust decomposition tool, as a 
combination of naked stellar photospheres and Planck curves of 
different temperatures, instead of the adapted Engelke function 
used by Sloan & Price ( 1998 ), we find a similar relation between 
the flux ratios. 



Fig. 11. The ZrO and TiO indices for the stars with emission 
features due to hydrocarbons, SiS or CO2. The other stars in the 
sample are shown as dots for comparison. 



The majority of the S-type stars in our sample have clas- 
sification SE1-SE4, and only a third of the stars have classifi- 
cation SE5-SE8, similar to the findings in Sloan & Price ( 1995 
1998 ). Based on our sample, we find that most stars from Group I 
can be classified as SE5-SE8 stars, while stars from Group II 
and Group III make up classes SE1-SE4. Although Group II 
and Group III stars clearly have different IR spectra at longer 
wavelengths, they have similar 10-12/mi features and hence 
there is no difference between the Sloan & Price classification 
of both groups. Hence, the stars with C/O ratios close to unity 
(Group III) make up the same IR classes as stars with low C/O 
ratios (Group II). This is contradictory to the idea that the 10- 
12 /mi emission is determined by the photospheric C/O ratio for 
AGB stars. However, there is still a likely difference between the 
position of these stars in Fig. 



uated below the powerlaw, w 



10 Most stars of Group III are sit- 



/hile the stars of Group II are above 
it. If this effect is real, we can conclude that the C/O ratio does 
not determine the Sloan & Price classification, but might explain 
the spread on this correlation, while the Sloan & Price class is 
determined by the mass-loss rate. 



7. Molecular emission features 

7.1. Hydrocarbon emission 

Hydrocarbon molecules exhibit very characteristic emission fea- 
tures in the 5-14//m region, the most prominent at 6.2, 7.9, 8.6, 
and 1 1 .2 /mi. These features arise from the bending and stretch- 
ing of the carbon and hydrogen bonds in the large molecules. 
Peet ers et al.| ( |2002) > showed that the 7.9 /mi PAH feature can 
be very different from source to source. They classified the PAH 
spectra in three classes, based on the peak wavelength of this 
feature. For type C sources, the peak of this feature is shifted to 
8.2 /mi, marking a clear difference with the 7.9 /mi feature of 
class B PAHs and the 7.6 /mi feature of class A PAHs. 



Smolders et al. ( 2010| l present the detection and identifica- 



tion of polycyclic aromatic hydrocarbon (PAH) emission in four 
S-type AGB stars from the present sample (BZ CMa, KR Cam, 
CSS 173 and CSS 757). Based on the peak-wavelength, BZ CMa 
is the only one that shows the B class PAHs resembling PNe, 
Herbig Ae/Be stars, and red supergiants ( |Tiele ns 2008). The 
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other three stars can be classified as rare class C sources resem- 



bling some of the PAHs found in post-AGB stars (Gielen et al. 
2009), the carbon-rich red giant HD 100764 (Sloan et al. 2007) 



and in a disk around the oxygen-rich K giant HD 233517 ( |Jura| 
|et al.|2006| ). 

Furthermore, these relatively red emission features are con- 
sistent with the strong correlation found between the centroid 
wavelength of the 7.9 /an feature and the temperature of the cen- 
tral star ( [Sloan et al.|2007) |Keller et al.|2008) |Acke et al.|20T0) ). 
The detection of PAHs in four S-type stars extends this corre- 
lation towards lower temperatures and more redshifted features. 
This is consistent with the hypothesis that Class C PAH spec- 
tra arise from a mixture of aliphatic and aromatic hydrocarbons, 
perhaps as isolated PAH molecules embedded in a matrix of 
aliphatic hydrocarbon bonds, found around stars with weak UV 
radiation fields (see |Pendleton & Alla mandola 2002 Fig. 17). 
The hydrocarbons around CSS 757, KR Cam, and CSS 173 thus 
represent the composition as condensed in the AGB wind, be- 
fore entering the interstellar medium where harsh UV radiation 
alters their chemical structure. 



In Fig. 1 1 we again show the ZrO and TiO indices of all stars 
in our sample. From this figure, it is clear that the stars with hy- 
drocarbon emission have C/O ratios that are higher than average. 
This is consistent with the idea that in stars with a nearly equal 
amount of carbon and oxygen atoms, non-equilibrium, shock- 
driven chemistry can form C2H2 molecules in the outer regions 
of the stellar outflow. These molecules are, further out in the 
wind, the necessary building blocks for the formation of hydro- 
carbons (Allamandol a et al.|1989| ). 

7.2. SiS emission 



Sloan et aL] ( |201 1[ ) identify the emission features at 6.6 and/or at 
13.5 /mi as molecular emission of SiS. These features appear at 
precisely the wavelengths where they are detected in absorption 
(Cami et al. 2009 ) . The lack of any wavelength shift strongly 
suggests that the emission is from SiS in the gas phase. Sloan 
et al. ( 20 1 1 ) l only discuss the two stars in our sample with the 
strongest emission features (CM Cyg and CSS 1005), an addi- 
tional five stars show a significant emission band at either 6.6 
or 13.5 /mi, all of which are classified as Group III stars (see 
Figs. [T2] ). The presence of these features in emission show that 
the SiS molecules are not only abundant in the stellar atmo- 
sphere, but also in the outer layers of the AGB wind. 

In Sect. [5] and 6.3.2 we have shown that both magnesium 
sulfide grains and SiS molecules are only observed around S- 
type stars with C/O ratios close to unity (see Fig. |4] and Fig. [9]). 
Therefore, this strong relation between the presence of SiS emis- 
sion and the dust classification could be explained as a shared 
dependence on the C/O ratio. However, it is possible that the 
magnesium sulfide grains in these stars are formed by a reaction 
of SiS molecules with Mg, as discussed in Sect. 8.4 



7.3. C0 2 emission 

The narrow features at 13.9, 15 and/or 16.2 /mi, visible in 9 ob- 
jects have been identified as CO2 emission lines (more specifi- 
cally 12 C 16 02). These lines are observed in approximately 30% 
of oxygen-rich AGB stars, most of which also show a 13 //m fea- 
ture Pusttanont et al.|1998[ [Sloan et al.|2003l . 

To extract these emission features, we fit a line under the 
features (at 13.7 and 14.1 /mi for the 13.9 /mi feature, 14.7 and 
15.2//m for the 15 /mi feature, and 15.9 and 16.5/mi for the 




15 20 25 

Wavelength [;im] 



Fig. 12. The spectra of all S-type stars in our sample with sus- 
pected SiS emission. The grey areas indicate the position of the 
6.6 and 13.5 SiS bands. At the bottom, CSS 100, one of the stars 
with SiS absorption bands, is shown for comparison. 



16.2 /mi feature). As a conservative criterion, we consider a de- 
tection of CO2 emission if all three features are significant at a 
99% confidence level. In our sample of S-type stars, we detect 
significant CO2 emission features in the 5 objects shown in Fig. 
13 The grey bands indicate the location of the emission features. 
It is clear from Fig. 
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that all stars with a positive detec- 
tion of CO2 emission lines have low C/O ratios. Because CO2 
requires an oxygen-rich environment to form, this dependence 
on the C/O ratio is not unexpected. Furthermore, although the 
signal-to-noise ratio and the resolution of the spectroscopic data 
do not allow for a quantitative analysis of the relationship be- 
tween the strength of the CO2 features and the 13 /mi feature, 
we can confirm that the 13 /mi feature is present in most of the 
stars with CO2 emission (4 out of 5 objects). This relation be- 
tween the presence of the 13 /mi emission and the detection of 
CO2 could be explained as a shared dependence on the C/O ratio. 

8. Discussion 

8. 1 . Carbon- versus oxygen-chemistry 

We find that almost all S-type AGB stars in our sample show 
either molecular absorption of molecules such as H2O or SiO, 
and/or emission due to dust species such as amorphous alumina 
or silicates. This mainly oxygen-dominated chemistry in the S- 
type stars is also confirmed by the stellar parameters, derived 
from the marcs model atmospheres, where we find that marcs 
models with C/O ratios between 0.5 and 0.99 can explain the 
optical spectra of all stars in our sample. This lack of carbon- 
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Fig. 13. The residual flux, after subtracting a line, of all stars 
with a positive detection of the CO2 emission features at 13.9, 
15 and/or 16.2 jum, indicated here with the grey bands. 



rich S-type stars is a direct consequence of the spectral classi- 
fication: for an AGB star to be classified as an S-type star, the 
ZrO bands need to be present in the optical spectrum. Ferrarotti 
& Gail ( 2002) show that the presence of ZrO in the stellar at- 
mosphere critically depends on the C/O ratio of the star. When 
the C/O ratio becomes larger than unity, the abundance of ZrO 
in the stellar atmosphere drops several orders of magnitude and 
hence the strength of the absorption bands drops, until the ZrO 
bands are not detectable anymore and the star is not classified as 
an S-type star ( |Piccirillo|1980 1. 

However, some stars additionally show absorption bands of 
SiS, hydrocarbon emission features or the broad emission fea- 
ture of magnesium sulfide dust, which are typical for carbon-rich 
sources. However, there are no stars in our sample with an IR 
spectrum that only shows carbon-rich characteristics. Because 
all stars in our sample have oxygen-dominated stellar atmo- 
spheres, it is remarkable that we can observe emission fea- 
tures that are generally only found in carbon-rich environments. 
The presence of these emission features indicates the (possibly 
strong) importance of non-equilibrium chemistry in the circum- 
stellar environment of these AGB stars. The best example of this 
is the presence of hydrocarbons around KR Cam, CSS 173, CSS 
757, and BZ CMa. The thermal equilibrium calculations pre- 
sented in Ferrarotti & Gail ( 2002) do not show the production 
of hydrocarbon molecules, because the majority of the carbon 
atoms would be locked in the stable and abundant CO molecule. 
Two main scenarios have been proposed to explain the pres- 
ence of the hydrocarbon emission in these S-type stars. The first 
scenario assumes the formation of hydrocarbons from the basic 



atoms, in the AGB wind itself. |Cherchneff| ( |2006| ) shows that, 
in AGB stars with a C/O ratio close to unity, non-equilibrium 
chemistry predicts the formation of C2H2 molecules in the stel- 
lar outflow ( |Cherclmeff eTaT1[T992l |Cau||2002|i, which are th e 
building blocks for hydrocarbons (Allamandola et al. 1989] >. 
Although this mechanism can explain the formation of hydro- 
carbons around S-type stars, C2H2 has a very distinct and sharp 
absorption feature at 13.7 /jm, which is not observed in our S- 
type stars. The second scenario starts from larger, carbonaceous 
grains that are broken down into smaller structures, such as hy- 
drocarbons. However, this scenario still requires the presence of 
carbonaceous dust grains in an oxygen-rich stellar environment. 



8.2. Variability characteristics versus dust properties 

The stellar winds in AGB stars are commonly attributed to the 
coupled effect between large amplitude pulsations and radiation 
pressure on dust grains ( Hoefner et al.|1998) . This relationship 
between the pulsations of the stellar atmosphere and the cir- 
cumstellar environment can be translated into in a relationship 
between the properties of the infrared spectrum and the vari- 
ability characteristics of the star. For oxygen-rich stars, a sur- 
vey of the globular cluster 47 Tuc has shown that not only the 
mass-loss rate, but also the dust mineralogy changes as the stars 
evolve towards higher luminosities and longer pulsation periods 
(|Lebzelter et al.|[2006l |van Loon et~aLl|20U6l |McDonald etaL] 
12011 p . 

We were able to derive the period and amplitude of the pul- 
sations for 55 out of the 87 S-type stars in our sample. This gives 
us a small subset of stars with and without a dusty circumstellar 
environment with known variability characteristics. An analysis 
of this subset shows no significant difference between the peri- 
ods and amplitudes of the stars with and without dust emission. 
Apparently for these S-type AGB stars, the period and ampli- 
tude of the pulsations are not enough to explain why approxi- 
mately 65% of the stars have no circumstellar dust, while the 
other 35% do. Furthermore, it is remarkable that we have de- 
tected 5 Mira-like pulsators without significant dust emission 
features. Looking at the slope of the IR spectra, we can argue 
that three of these stars have a continuum excess, which might 
arise from amorphous carbon or metallic iron (AU Car, RX Car 
and CSS 724). However, CSS 794 and XY CMi show no dust 
excess at all. 

Furthermore, a multivariate analysis of the variability char- 
acteristics, the stellar parameters and the results from the dust 
decomposition tool found no significant correlation with either 
amplitude or period. From this, we can conclude that during this 
short evolutionary phase where a star is visible as an S-type star, 
the period and amplitude are not indicative for the mineralogy or 
the dust mass-loss rate. 



8.3. Dust formation around S-type AGB stars 

In Sect. [6] the dust decomposition tool was used to quantify the 
contribution of dust emission features to the infrared flux be- 
tween 5 and 35 fim. A multivariate analysis of these parame- 
ters, together with the stellar parameters (r e ff, C/O and [s/Fe]) 
yields four significant correlations (using a 99% certainty in- 



terval) shown in Fig. 15 As an indication for significance of 



the correlation, we give the probability of obtaining a Pearson, 
Spearman or Kendall correlation coefficient at least as extreme 
as the one that was actually observed, assuming that there is no 
correlation between the two given variables. In order to interpret 
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Fig. 14. Schematical overview of the proposed idea for dust for- 
mation in S-type AGB stars. 

these p-values correctly, one has to keep in mind that Pearson 
assumes a linear correlation, while Spearman and Kendall are 
ranking coefficients. Because we look at ten parameters and thus 
effectively investigate 45 possible correlations, we have to ac- 
knowledge the possibility for a false-positive correlation. All im- 
portant correlations are shown in Fig. [15] 

Silicates and total dust emission If we remove the outlier 
V376 Aur, we find a significant correlation between the to- 
tal dust emission and the flux percentage of silicate fea- 
tures, shown in Fig. 15(a) This correlation indicates that 
stars with a larger dust excess, and thus a higher dust mass- 
loss rate, also tend to form silicates in their circumstellar 
environment. This is in close agreement with the study of 
Galactic bulge AGB stars ( |Blommaert et al. |2006 1 and the 
study of AGB stars in the Large Magellanic Cloud presented 
i Dijkstra et al.| ( |2005) . Here, the authors find that the dust 



features displayed by oxygen-rich stars are mainly ascribed 
to alumina for low dust mass-loss rates, while stars with 
higher dust mass-loss rates show an increasing amount of 
silicates. Furthermore, if we look at the total dust excess in 
the three groups, we find that, on average, Group III stars 
have a higher dust excess than Group I and Group II stars. 
Although the difference is only tentative, this is still a note- 
worthy dichotomy between the three groups, because this is 
what could be expected based on the results of |Blommaert| 



|etaL]p006) l and |Dijkstra et al.|(|2005| >. 
Silicates and alumina In Fig. 15(b) we show the significant 
and strong inverse correlation between the relative flux con- 
tribution of alumina and silicates. This correlation can be 
explained by the standard dust condensation sequence, if we 
take into account freeze-out. Alumina seeds are formed early 
on in the dust-condensation sequence at high temperatures, 
but if the density in the wind is high enough, eventually sil- 
icates will form and hence the relative flux contribution of 
alumina will decrease as the flux contribution of silicates in- 
creases. 

Magnesium sulfides and alumina In Fig. 



15(c) a strong cor- 



relation between the relative flux contribution of magnesium 
sulfides and alumina is shown. This correlation is unexpected 
because the presence of alumina grains points to an oxygen- 
rich circumstellar chemistry while magnesium sulfide emis- 
sion is usually found in carbon-rich stars. Because this is the 



first time that this correlation is observed, it can be consid- 
ered as a new constraint on any dust formation hypothesis 
for S-type AGB stars. 
Magnesium sulfides and silicates In Fig. 



15(d) we show the 



inverse correlation between the relative flux contribution of 
silicates and magnesium sulfides. This correlation does not 
add extra constraints, because it is a consequence of the cor- 
relations between the relative flux contribution of (i) silicates 
and alumina and (ii) magnesium sulfides and alumina. 

We propose that the dust formation in S-type stars is crit- 
ically dependent on the C/O ratio and the total mass-loss rate. 
Because we expect that silicates only form further out in the 
wind where the temperature is low enough, the dust condensa- 
tion sequence can be cut short before silicates start to condense. 
This could happen because (i) the density in the wind is too low 
to initiate the necessary chemical reactions (freeze-out) or (ii) 
there is not enough free oxygen to form the next dust species 
in the dust condensation sequence (oxygen depletion). At this 
point, stars with C/O ratios close to unity can only form alu- 
mina grains and possibly magnesium sulfide. Stars with lower 
C/O ratios still have enough free oxygen atoms to form alumina 
and the carriers of the 13 fim and 19.8 /urn features. In stars with 
higher mass-loss rates, the density throughout the wind is high 
enough to start forming silicates. Because the formation of sil- 
icates requires a large number of oxygen atoms, the density at 
which silicates can condense might depend on the C/O ratio. 

This hypothesis can explain perfectly why we can find 
Group I stars with any C/O ratio, as seen in Fig.[9| while Group II 
stars are clearly oxygen-rich and Group III stars have C/O ratios 
close to unity. Furthermore, this hypothesis explains the differ- 
ence in total dust excess between the groups and it takes into 
account the tentative correlation between the strength of the sil- 
icate emission and the total dust excess. 



8.4. The formation of magnesium sulfides in S-type AGB 
stars 

A number of observational and theoretical studies on the 26- 
30 /im emission feature in AGB and post- AGB stars, and in the 
environment of planetary nebulae show that magnesium sul- 
fide is a wide-spread dust component in carbon-rich evolved 
stars (|Hony et al.||2002"l |Volk et al.||2"002"l |Zijlstra et ai^|2006] 
|Lagadec et al.|2007) . As magnesium sulfide is usually observed 
in circumstellar environments dominated by carbon chemistry, 
the chemical models have focused on magnesium sulfide for- 
mation in a carbon-rich environment with and without silicon 



carbide condensation (Zhukovska & Gail 2008). There are two 
main pathways to form magnesium sulfide grains. If the sulfur is 
locked in SiS molecules, the necessary reaction is: 

Reaction I : SiS + Mg -> MgS (solid) + Si 

If the sulfur is free from the SiS molecules because, for ex- 
ample, the silicon is consumed in the formation of silicon car- 
bide, the free sulfur forms H2S molecules and the main reaction 



Reaction II : H 2 S + Mg -> MgS (solid) + H 2 



Zhuk ovska & Gail| (2008) find for carbon-rich AGB stars that 
only the growth of magnesium sulfide as mantles on SiC cores, 
following reaction II, can produce enough magnesium sulfide 
dust to explain the observations. Hence, following this forma- 
tion mechanism, the production of magnesium sulfide depends 
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Fig. 15. All significant correlations presented in Section 8.3 Above each figure, you can find the Spearman, Kendall and Pearson 



probabilities of each of the correlations. Please keep in mind that Pearson can only be used for linear relationships, while Spearman 



and Kendall are both rank-test and thus can be used for any relation between two parameters. The correlations in figures (b) (c) and 



(d) are clearly significant, while the correlations in figure (a) is only tentative 



on silicon carbide grains (i) to free the sulfur from SiS molecules 
(so that H2S can be formed) and (ii) as core grains on which the 
magnesium sulfide can form a mantle. Since we do not see strong 
silicon carbide emission features in any of our stars (although we 
cannot completely rule out the presence of a weak emission fea- 
ture in the 10-12 /jm range) we conclude that alternative seeds 
are necessary. 

If reaction II is the dominant formation path to form mag- 
nesium sulfides, we need to free the sulfur from SiS before this 
reaction can take place. In a mixed chemistry region, there is a 
much wider variety of silicon-rich dust species as compared to a 
purely carbon-rich chemistry (see Fig. [5]). Already at relatively 
high temperatures the corundum reacts with circumstellar gas 
to form gehlenite and further interactions lead to other silicon- 
bearing dust grains such as akermanite, diopside, anorthite, and 
eventually silicate grains. All these dust species can consume sil- 



icon, hence allowing the free sulfur to form the H2S molecules, 
necessary for reaction II. 

Although these silicon-bearing grains have emission peaks 
in the 10-20/vm range, the infrared spectra of most of the stars 
with clear magnesium sulfides emission also show amorphous 
alumina, without clear observational evidence for any silicon- 
bearing dust species. Hence, the observations indicate that the 
presence of silicates or other silicon-bearing dust species is not 
a strict requirement for the formation of magnesium sulfides. 

Because many stars with magnesium sulfide emission show 
additional emission peaks at 6.6 and 13.5 fim, due to molecular 
emission of SiS, the observations seem to indicate that reaction 
I is the main formation mechanism in these S-type AGB stars. 
Because sulfur is less abundant than magnesium, we can expect 
that more influx of SiS molecules in the wind would lead to 
more magnesium sulfides being formed. Furthermore, because 
SiS molecules are only formed in S-type stars with C/O ratios 
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larger than 0.96, it automatically follows that we will only detect 
magnesium sulfide dust in stars with similar high C/O ratios, as 

is shown in Fig. [9] 

According to Zhukovsk a & Gaii] ( |2008| l, the formation path 
of magnesium sulfide is more efficient as heterogeneous growth 
on a different type of dust particles that already have moderate 
radii. The dust species that may serve as carrier grains need to 
form dust grains with a size not much smaller than about 0. 1 yum 
and need to be formed early on in the wind. The correlation 
shown in Fig. 



We also find a strong relation between the C/O ratio and 
mineralogy. We propose that, for stars with low mass-loss rates, 
the dust-condensation sequence is cut short and the remain- 
ing amount of free oxygen, together with the presence of SiS 
molecules in the wind determines whether sulfides will be 
formed (Group III) or whether we will observe a 13yimi and 
20 /im feature (Group II). Stars with higher mass-loss rates 
can complete the dust formation sequence and form silicate 
grains (Group I). This scenario is consistent with the correlations 



15(c)|indicates an (indirect) relation between alu- shown in Fig. 15 



mina grains and the magnesium sulfide formation. Because alu- 
mina can form very early on in the wind (at temperatures above 
1700 K) and because of the correlation between the relative flux 



contribution of alumina and magnesium sulfides (see Sect. 8.3 1 



we propose the hypothesis that magnesium sulfides in S-type 
AGB stars are formed on amorphous alumina grains. In addi- 
tion, this formation mechanism can explain the presence of SiS 
emission bands at 6.6 and/or 13.5 /jm in many Group III stars, 
because it is a key ingredient for the formation of magnesium 
sulfides. 



9. Results and conclusions 

We present a new sample of Spitzer-IRS spectra for 87 galactic 
S-type AGB stars which were selected based on the presence of 
a ZrO absorption band at optical wavelengths, pulsation charac- 
teristics and K-band magnitude. Although the sample is slightly 
biased towards lower mass-loss rates, we show that we exclude 
less than 5% of the stars with high mass-loss rates. In addition 
to the IR spectra, we obtained optical spectra (low and/or high 
resolution) and infrared photometry for the stars in our sample. 

We derived the chemical parameters (T e ff, C/O, [s/Fe] and 
[Fe/H]) for a large subset of the sample, by comparing the avail- 
able optical spectra to marcs model atmospheres. 

For a large subset of the sample, we provide pulsation char- 
acteristics, based on the ASAS, NS VS and AAVSO light curves. 
The sample contains 6 irregular pulsators, 38 semi-regular vari- 
ables and 21 Miras. It is remarkable that 5 of the 17 Miras do 
not show any dust emission features. This does not imply that 
these stars have no dusty wind, because we cannot exclude the 
presence of dust species, such as metallic iron and amorphous 
carbon. 

We find that the presence of the SiS absorption band in the S- 
type stars without dust emission features is critically dependent 
on the C/O ratio. Stars with C/O ratios lower than 0.96 do not 
show significant absorption at 6.6 or 13.5 jum while all naked 
stellar atmospheres with C/O ratios above 0.96 show significant 
absorption bands. This opens up the possibility to use the SiS 
absorption band as an indicator for the C/O ratio. 

We present a positive detection and identification of PAH 
emission in 4 S-type AGB stars, a detailed description and dis- 
cussion is given by Smolders et al. ( 2010| l. 

We use a simple tool to decompose the IR spectra of the dust- 
rich stars. This tool is used to identify the dust species that are 
visible in the IR spectra and to quantify the relative flux con- 
tribution of each dust type to the total flux. Based on this dust 
decomposition tool, we are able to classify the IR spectra into 
three distinct groups: (i) stars with strong silicate/gehlenite fea- 
tures, (ii) stars with strong 13/im and 20 //m emission features 
and (iii) stars with magnesium sulfides and amorphous alumina. 
The fact that we find three separate classes is already a remark- 
able result, for example, no star has both a 13yum feature and 
magnesium sulfide emission. 



Finally, we propose the hypothesis that magnesium sulfides 
in S-type stars are formed by a reaction of SiS molecules with 
Mg. This is based (i) on the observation that magnesium sulfides 
are only observed in stars with C/O ratios close to one, where 
SiS molecules are formed and (ii) on the detection of SiS emis- 
sion in half of the stars with magnesium sulfide emission, indi- 
cating that the SiS molecules are present in the circumstellar en- 
vironment. Furthermore, based on the formation mechanism of 
magnesium sulfides in carbon-rich AGB stars, we propose that 
magnesium sulfide is formed as mantles around dust grains, for 
example amorphous alumina. 
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Table 3. Effective temperature, C/O ratio and gross s-process overabundance of S stars. Hen 4-NNN lists S stars as in the Simbad 
database from an unpublished list of Henize (see e.g. Van Eck et al.|2 000 1.The "uncertain" label in the comment column means that 
the convergence to the quoted set of physical parameters was not as consistent as in all the other cases. 





left 




C/O 


[s/Fe] 


1 nmmpnte 
I^OIIIIIICIIIS 


127 


3300 ± 


100 


0.95 


+ 


0.025 


1.5 


± 


0.5 




166 


3300 ± 


100 


0.98 


+ 


0.010 


1.7 


+ 


0.3 




173 


3600 ± 


100 


0.97 


± 


0.020 


1.5 


± 


0.5 


uncertain 


177 


3350 ± 


150 


0.96 


± 


0.010 


1.0 


± 


0.5 


KR Cam 


179 


3250 ± 


50 


0.98 


± 


0.010 


1.5 


+ 


0.5 




191 


3200 ± 


100 


0.84 


+ 


0.088 


1.5 


± 


0.5 


GH Aur ; H„ emission ; Mira ; uncertain 


196 


3250 + 


50 


0.84 


± 


0.088 


0.5 


+ 


0.5 


V376 Aur 


202 


3500 ± 


100 


0.70 


± 


0.200 


1.0 


± 


0.3 


CD-39 2449 ; Hen 4-8 ; 


220 


3300 ± 


100 


0.98 


+ 


0.010 


1.5 


± 


0.5 


CF Gem 


267 


3350 ± 


50 


0.99 


+ 


0.010 


1.5 


+ 


0.5 




281 


3300 ± 


100 


0.99 


± 


0.010 


1.0 


+ 


0.5 


BZ CMa 


282 


3300 ± 


100 


0.97 


± 


0.020 


1.5 


+ 


0.5 




315 


3450 + 


50 


0.63 


+ 


0.125 


0.6 


+ 


0.4 


CD-43 2976 


325 


3200 ± 


100 


0.71 


+ 


0.210 


0.5 


± 


0.5 


XY CMi 


328 


3450 ± 


150 


0.93 


± 


0.025 


1.0 


± 


0.3 


CD-30 4223 ; Hen 4-19 ; 


352 


3400 ± 


100 


0.70 


± 


0.200 


1.3 


± 


0.3 




361 


3300 ± 


100 


0.97 


+ 


0.020 


1.5 


± 


0.5 


Hen 4-23 


371 


3400 ± 


100 


0.98 


+ 


0.010 


1.5 


+ 


0.5 


Hen 4-26 ; uncertain 


380 


3150 ± 


250 


0.98 


± 


0.010 


1.0 


+ 


0.7 


Hen 4-27 ; H B emission 


393 


3450 + 


50 


0.97 


± 


0.020 


1.5 


+ 


0.5 


Hen 4-29 


417 


3450 + 


50 


0.98 


+ 


0.010 


1.5 


+ 


0.5 


H tt emission 


427 


3200 ± 


100 


0.97 


+ 


0.020 


1.5 


± 


0.5 




431 


3300 ± 


100 


0.96 


± 


0.010 


1.5 


+ 


0.5 


AO Gem ; H ff emission 


438 


3150 ± 


150 


0.96 


± 


0.010 


1.5 


+ 


0.5 




443 


2900 ± 


200 


0.71 


+ 


0.210 


0.5 


± 


0.5 


EX Pup ; H ff emission 


453 


3400 ± 


100 


0.70 


+ 


0.200 


0.5 


+ 


0.5 




454 


3400 ± 


100 


0.71 


± 


0.210 


1.0 


± 


0.2 




457 


3450 ± 


50 


0.63 


± 


0.125 


0.5 


+ 


0.5 


Hen 4-38 ; uncertain 


472 


3350 + 


50 


0.63 


+ 


0.125 


0.1 


+ 


0.1 


MS star 


478 


3400 ± 


100 


0.71 


+ 


0.210 


0.5 


± 


0.5 




480 


3350 + 


50 


0.63 


± 


0.125 


0.4 


± 


0.3 




489 


3450 + 


50 


0.63 


± 


0.125 


1.0 


± 


0.3 


CD-29 5912 ; Hen 4-44 


528 


3450 + 


50 


0.83 


+ 


0.075 


0.4 


± 


0.4 




531 


3250 + 


50 


0.91 


+ 


0.012 


1.2 


+ 


0.2 




553 


3300 ± 


100 


0.70 


± 


0.200 


0.4 


+ 


0.3 




558 


3300 ± 


100 


0.71 


± 


0.213 


0.9 


± 


0.3 


H„ emission 


584 


3300 ± 


100 


0.99 


+ 


0.010 


1.5 


+ 


0.5 


Hen 4-62 ; uncertain 


598 


3250 ± 


150 


0.97 


+ 


0.010 


1.5 


± 


0.5 




603 


3400 ± 


100 


0.70 


± 


0.200 


0.4 


+ 


0.4 


CD-30 7391 


604 


3300 ± 


100 


0.98 


± 


0.010 


1.5 


+ 


0.5 




613 


3450 + 


50 


0.98 


+ 


0.010 


1.7 


± 


0.3 




622 


3200 ± 


100 


0.99 


+ 


0.010 


0.5 


+ 


0.4 




633 


3300 + 


50 


0.91 


+ 


0.012 


0.2 


± 


0.1 


Hen 4-73 ; uncertain 


635 


3400 ± 


100 


0.98 


± 


0.010 


1.5 


+ 


0.5 


H a emission 


640 


3250 + 


50 


0.97 


± 


0.010 


1.0 


+ 


0.3 


Hen 4-75 


650 


3450 + 


50 


0.70 


+ 


0.200 


0.8 


± 


0.5 


Hen 4-77 


657 


3300 ± 


100 


0.97 


+ 


0.010 


1.5 


+ 


0.5 


Hen 4-81 ; H„ emission 


665 


3300 ± 


100 


0.97 


± 


0.020 


1.5 


+ 


0.5 




682 


3400 ± 


100 


0.71 


± 


0.210 


0.8 


+ 


0.3 


Hen 4-92 


690 


3150 ± 


150 


0.71 


+ 


0.210 


0.3 


+ 


0.3 


RX Car ; H ff emission ; Mira ; uncertain 


718 


3300 ± 


100 


0.71 


+ 


0.213 


0.5 


± 


0.5 


Hen 4-106 ; uncertain 


719 


3400 ± 


100 


0.99 


± 


0.010 


1.5 


+ 


0.5 




739 


3200 ± 


100 


\J.yl 


+ 


U.Ulz 


n i 
U.l 


+ 


U. 1 


Mira 


757 


3300 ± 


100 


0.98 


+ 


0.010 


1.5 


± 


0.5 




763 


3300 ± 


100 


0.71 


+ 


0.213 


0.5 


+ 


0.5 




783 


3450 ± 


150 


0.71 


± 


0.210 


0.5 


+ 


0.5 




794 


3150 + 


150 


0.96 


± 


0.010 


0.5 


+ 


0.5 


Hen 4-122 


806 


3650 + 


50 


0.83 


+ 


0.075 


0.5 


+ 


0.5 


high extinction; uncertain 


881 


3250 ± 


150 


0.99 


+ 


0.010 


1.0 


± 


0.5 




921 


3350 + 


50 


0.63 


± 


0.125 


0.3 


± 


0.3 




1005 


3350 ± 


150 


0.97 


± 


0.020 


1.5 


± 


0.5 




1180 


3300 ± 


100 


0.99 


+ 


0.010 


1.0 


± 


0.3 


CM Cyg ; strong H ff emission 


1305 


3250 ± 


150 


0.70 


+ 


0.200 


0.3 


± 


0.3 


CE Lac 


1336 


3300 ± 


100 


0.63 


± 


0.125 


0.4 


+ 


0.3 
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Appendix A: Results of the dust decomposition 
tool 
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CSS 1: x 2 =48.18 
Silicate: 3.82±4.31 
Gehlenite: 1.21 ±1.50 
Alumina: fi8.96±6.0G 
Spinel: 0.02 ±0.10 
Mgl0Fe90O: 2.84±2.32 
Mg90FelOS: 23.15±6.07 
Total fit 




RX Psc: x 2 =14.04 
Silicate: 22.63±8.71 
Gehlenite: 52.05±3.67 
Alumina: 19.57 ±6,12 
Spinel: 0.00±0.00 
Mgl0Fe90O: 0.00±0.00 
Mg90FelOS: 5.75±2.52 
Total fit 



15 20 25 

Wavelength [micron] 




Wavelength [micron] 



NX Per: x 2 -18.75 
Silicate: 20.81±3.30 
Gehlenite: 19.49±l.f>3 
Alumina: 51.01±3.22 
Spinel: 0.62±0.51 
Mgl0Fe90O: 2.21±1.23 
Mg90FelOS: 5-86 ±2.21 
Total fit 




GH Aur: x =13.20 
Silicate: 4.50 ±4.44 
Gehlenite: 2.42±2.17 
Alumina: 57.10 ±7.18 
Spinel: 4.34±1.14 
Mgl0Fe90O: 14.06±2-24 
Mg90FelOS: 7.57±3.04 
Total fit 




Wavelength [micron] 



Wavelength [micron] 



V376 Aur: x 2 -82.54 
Silicate: 20.94±3.28 
Gehlenite: 16.82±1.61 
Alumina: 43.90 ±2.13 
Spinel: 2.34 ±0.52 
Mgl0Fe90O: 11.42 ±1.06 
Mg90FelOS: 4.58 ±1.5.0 
Total fit 




CD-39 2449: x 2 -59.17 
Silicate: 83.71 ±2.53 
Gehlenite: 13.41 ±1.39 
Alumina: 0.06±0.41 
Spinel: 0.00±0.00 
MglOFe900: 1.71±0.80 
Mg90FelOS: 1.10±1.09 
Total fit 




Wavelength [micron] 



Wavelength [micron] 



Fig.A.l. The results of the dust decomposition tool for all stars with dust features. Top panel: the black line shows the Spitzer 
spectrum, the red dotted line the fit the underlying estimate for the stellar continuum, as derived from the dust decomposition tool 
and the solid, dark blue line indicates the total fit to the spectrum. Solid lines shown in the legend indicate significant contributions 
to the spectrum, the dotted lines the non-significant contributions. Bottom panel: the residuals after subtracting the total fit. 
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HV 8057: X 2 -15.52 
Silicate: 22.79±4.62 
Gehlenite: 20.87±2.17 
Alumina: 50.62±3.68 
Spinel: 0.50±0.57 
Mgl0Fe90O: 2.28±1.10 
Mg90FelOS: 2.93±2.14 
Total fit 




CSS 361: x 2 =16.55 
Silicate: 18.82±5.92 
Gehlenite: 38.25±2.53 
Alumina: 37.02±4.37 
Spinel: 2.38±0.91 
Mgl0Fe90O: 0.00±0.01 
Mg90FelOS: 3.53±1.10 
Total fit 
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CSS 380: x 2 =24.03 
Silicate: 10.38±5.94 
Gehlenite: 10.58±2.46 
Alumina: 66.33±3.26 
Spinel: 0.00 ±0.01 
Mgl0Fe90O: 0.17±0.43 
Mg90FelOS: 12.54±2.79 
Total fit 




CSS 427: x 2 =51.40 
Silicate: 23.96±4.73 
Gehlenite: 12.28±2.36 
Alumina: 45.03±4.96 
Spinel: 0.09 ±0.25 
Mgl0Fe90O: 11.06 ±4.13 
Mg90FelOS: 7.58±4.18 
Total fit 
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AO Gem: x 2 =38.10 
Silicate: 4.74 ±4.58 
Gehlenite: 24.70±2.51 
Alumina: (>1.44±3.54 
Spinel: 0.04±0.20 
Mgl0Fe90O: 0.18±0.45 
Mg90FelOS: 8.90±2.34 
Total fit 




CSS 438: x 2 "17.73 
Silicate: 16.23±2.66 
Gehlenite: 18.31±0.81 
Alumina: 54.96±2.03 
Spinel: 0.27±0.27 
Mgl0Fe90O: 0.99±0.77 
Mg90FelOS: 9.25 ±1,57 
Total fit 
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A.2. Continued 



K. Smolders et al.: The Spitzer Spectroscopic Survey of S-type Stars, Online Material p 4 




Fig.A.3. Continued 
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A.4. Continued 
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CSS 749: x 2 =12.74 
Silicate: 7.81 ±6.41 
Gehlenite: 12.89±2.66 
Alumina: 57.76±10.78 
Spinel: 1.50 ±0.93 
Mgl0Fe90O: 9.07±2.06 
Mg90FelOS: 10.97±3.54 
Total fit 




CSS 763: x 2 =12.81 
Silicate: 12.62±8.89 
Gehlenite: 24.74 ±3.95 
Alumina: 49.83±11.84 
Spinel: 1.35±1.22 
Mgl0Fe90O: 3.27±2.31 
Mg90FelOS: 8.19±2.90 
Total fit 
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CSS 987: x 2 =60.50 
Silicate: 0.26±1.04 
Gehlenite: 10.59 ±1.49 
Alumina: 66.66±2.88 
Spinel: 0.00±0.04 
Mgl0Fe90O: 2.34±1.81 
Mg90FelOS: 20.14±4.72 
Total fit 




CSS 1005: x 2 =15.28 
Silicate: 0.00±0.03 
Gehlenite: 23.08 ±2.01 
Alumina: 59.42±6.33 
Spinel: 0.96±1.56 
Mgl0Fe90O: 0.8f)±1.95 
Mg90FelOS: 15.08 ±0.25 
Total fit 




Wavelength [micron] 



Wavelength [micron] 




CSS 1160: x 2 -47.03 
Silicate: 9.15±5.00 
Gehlenite: 22.03±2.58 
Alumina: 58.04±3.28 
Spinel: 0.00±0.03 
MglOFe900: 1.71 ±1.45 
Mg90FelOS: 9.08±3.54 
Total fit 
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A.5. Continued 



K. Smolders et al.: The Spitzer Spectroscopic Survey of S-type Stars, Online Material p 7 





Wavelength [micron] 




A.6. Continued 



